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Abstract. Spectroscopic observations of the 2006 outburst of RS Oph at 
both infrared (IR) and X-ray wavelengths have shown that the blast wave has 
decelerated at a higher rate than predicted by the standard test-particle adia- 
batic shock-wave model. The observed blast-wave evolution can be explained, 
however, by the diffusive shock acceleration of particles at the forward shock 
and the subsequent escape of the highest energy ions from the acceleration re- 
gion. Nonlinear particle acceleration can also account for the difference of shock 
velocities deduced from the IR and X-ray data. We discuss the evolution of the 
nova remnant in the light of efficient particle acceleration at the blast wave. 



Introduction 



X-ray observation s of the latest outburst of RS Oph ( Sokoloski et al. 1 120061 ; 
iBode et al. 20061 ) have allowed to clearly identify the forward shock wave ex- 
panding into the red giant (RG) wind and to estimate the time evolution of its 
velocity, v s , through the well-known relation for a test-particle strong shock: 
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where k is the Boltzmann constant, T s is the measured postshock temperature 
and \xrrtii is the mean particle mass. The X-ray emission has revealed that after 
an ejecta-dominated, free expansion stage lasting ~6 days, the remnant rapidly 
evolved to display behavior characteristic of a shock experiencing significant 
radiative cooling. The lack or the very short duration of an adiabatic, Sedov- 
Taylor phase differs from the remnant evolut ion model developed after the 1985 
outburst (see lO'Brien. Bode, Kahn|[l9 92. and references therein). 

The time-dependence of shock velocity has also been m easured by IR spec- 
troscopy (|Das. Banerjee. fc Ashok |[2006l : lEvans et al. 1120071 ) . Although the gen- 
eral behavior of the shock evolution was found to be consistent with that deduced 
from the X-ray emission, the shock velocities determined from the IR data are 
significantly greater than those obtained using equation (1) together with the 
X-ray measurements of T s (see Fig. la). 

We have shown that production of nonthermal particles b y diffusive shock 

acceleration at the blast wave can account for these observations (Tatischeff fc Hernanz 
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2007). We have proposed that the difference of shock velocities deduced from the 



IR and X-ray data is due to the use of equation (1), which is known to underes- 



timate v s when particle ac celeration is efficient (jDecourchelle. Ellison. Sz Ballet 



l2000l : lEllison et al. 1 120071 ) . The observed early cooling of the blast wave is 
explained by the escape of high-energy particles from the acceleration region, 
which we found to be much more effective to cool the shock than radiative losses. 
Thus, nonlinear diffusive shock acceleration has important implications for the 
evolution of the nova remnant. 



2. The density in the red giant wind 

The RG wind density can be estimated from both the photoelectric absorption 
of the postshock X-rays and the free-free absorption of the radio synchrotron 
emission, which presumably arises from electron acceleration at the blast wave. 
Using the absorbing colu mn density Njj measured with the Swift X-Ray Tele- 
scope (jBode et al. 1 12006). and assuming the abundances in the RG wind to be 
solar (see Wallerstein et al., these proceedings), we obtain from equation (4) of 
Bode et al. I (|2006l ): M RG /u KG « 4 x 10 g cm , where Mrg an d ifRG are the 



RG mass-loss rate and wind terminal speed, respectively. It is noteworthy that 
the values of Nh are s till uncertain, bec ause they were determined by rough 
single-temperature fits ( Bode et al. 1120061 ). 



To calculate t he optical depth to radio free- free absorption, we use the 
formula derived by lO'Brien. Bode. Kahnl (|l992l ) , which allows for the effects 
of H recombination in the stellar wind. For Mrg/urg = 4 x 10 13 g cm -1 and 
a uniform wind temperature T\y in the range (1-3) xlO 4 K, we find the time 
of free-free optical depth unity at 6.03 GHz to be 20-30 days. This result is 
consistent with the broad peak of radio emission observed with the Very Large 
Array (VLA) and MERLIN after t = 20 days, but is inconsistent with the early 
rapid rise of radio emission detected at days 4 and 5 after outburst (Eyres et 
al., these proceedings). If this early emission is associated with the blast wave 
expanding at v s ~ 4300 km s^ 1 (Fig. la), it implies Mr g/^rg < 7xl0 12 g cm -1 . 
From GMRT observations at frequencies < 1.4 GHz, Kantharia et al. ( 20071 ) 
obtained Mrg/'Urg = 5 x 10 12 g cm -1 . However, the X-ray- and radio-based 
estimates of the RG wind density could be reconciled if the first early peak of 
radio emission is associated with a synchrotron jet expanding much faster than 
the blast wave (see O'Brien et al. and Rupen et al., these proceedings). 



3. Properties of the Cosmic-Ray Modified Shock 

The effects of particle acceleration at the blast wave depend on the strength of 
the magnetic field just ahead of the shock, Bq. Assuming the magnetic field 
in the RG wind to be in equipartition with the thermal energy density, we 
have Bq = ob^Mrg/cTV /ujicr^fimH) ' 5 , where ob is a factor accounting for a 
possible amplification of the magnetic field in the shock region and r s is the shock 
radius. Here and in the following, we adopt 0.3 = 1, -MrgA*rg = 2x 10 13 g cm -1 
and T\y = 3 x 10 K, which give Bq = 0.06 G at t = t%. With such values of 
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Figure 1. (a) Time-dependence of forward shock velocity as deduced from 
FWZI of IR emission lines (filled circles: Das et al. 2006, filled triangles: 
Evans et al. 2007) and X-ray measurements of the postshock temperature 
(open triangles: Sokoloski et al. 2006, open squares: Bode et al. 2006). The 
IR data can be modeled by (solid line) v s (t) — 4300(i/ii) Q " km s _1 , where 
t\ = 6 days and a v = (-0.5) for t < t\ (t > t\). (b) Calculated postshock 
temperature for two values of n- m j compared to the RXTE and Swift data. 



Bq, we have shown ( Tatischeff Sz Hernanz 20071 ) that accelerated electrons and 
protons can have achieved TeV energies in few days after outburst. 

To study the modification of the shock structure induced by the backre- 
action of the energetic i ons, we use the mode l of n onlinear diffusive shock ac- 
celeration developed by Berezhko &: Ellison I ( 19991 ). with the prescription of 
Blasi. Gabici. fc Vannoni I (12005) for the particle injection. Calculated temper- 



atures of the postshock gas are shown in Fig. lb for two values of the parameter 
r/[ n j , which is the fraction of postshock thermal protons inj ected into the accelera- 
tion p rocess. The other parameters of the model are as in iTatischeff fc Hernanz 
((2003) ■ We see that the temperatures measured with RXTE and Swift can be 
well reproduced with 7^=1. 5x 10~ 4 . For r/i n j=10 -5 , the test-particle approxima- 
tion applies and the standard relation between v s and T s (eq. [1]) overestimates 
the temperature. The model also provides the shock compression ratio and ac- 
celeration efficiency. We find that the energy carried off by high-energy particles 
escaping the shock region represents 10%-20% of the total energy flux at t > t\ . 



4. Calculated X-Ray Fluxes 

We compare in Fig. 2 calculated X-ray fluxes w ith the RXTE and Swift data of 
Sokoloski et al. I ( 20061 ) and iBode et al. I ( 20061 ). respectively. Using theoretical 



emissivities from the MEKAL plasma emission model and calculated values of 
T s and the postshock density from the nonlinear shock acceleration model, we 
determine the characteristic emission volume, V(t) = 4irrg(t)Ar s (t), required to 
reproduce the measured fluxes. We assume the relative thickness of the X-ray 
emitting shell to be of the form Ar s (t)/r s (t) = f r (r s (t) /r s (t{)) ar , where f r and 
a r are free parameters. The two curves shown in Fig. 2 are for a r = and 
The first value corresponds to an adiabatic expansion of the shock 
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Figure 2. Time-dependence of absorption-corrected X-ray flux in the 0.7- 
10 keV energy range. Using the MEKAL plasma emission model, corrections 
ranging f rom 0.76 to 0.90 were a pplied to the RXTE PCA fluxes, which are 
given by ISokoloski et~aT~l (|2006fl in the 0.5-20 keV energy range. Dashed 
curve: a r = Q; solid curve: a r — a v (see text). 



wave. It is clear that the observed flux evolution implies a shrinking of the 
relative shell thickness at t > ti, as expected from a cooling shock. The nor- 
malization of the calculated curves to the data gives f r = 6.8% x (D^ pc /1.6) 2 , 
where -Dkpc is the source distance in kpc. This result is consistent with hydro- 
dynamic models of supernova remna nt (SNR) evolution with efficient cosmic- 
ray acceleration ( Ellison et al. 20071 ). as well as with SNR observations (e.g. 



Warren et al. 1120051 ). We note that a value of Mrg/^rg < 7 x 10 12 g cm l , as 



needed if the early radio emission is associated with the blast wave (see § 2) 
would require a too large X-ray emission volume to reproduce the measured 
fluxes (Fx oc V (Mug / urg) 2 ) ■ This gives some support for the assumption that 
the early radio peak is due to a fast synchrotron jet. 
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